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ABSTRACT 

We present results on the spectroscopic study of the ionized gas in the high redshift 
radio galaxy USS0828+193 at z —2.57. Thanks to the high S/N of the emission lines 
in the Keck spectrum, we have been able to perform a detailed kinematic study by 
means of the spectral decomposition of the emission line profiles. This study reveals the 
existence of two types of material in this object: a) a low surface brightness component 
with apparent quiescent kinematics consistent with gravitational motions and b) a 
perturbed component with rather extreme kinematics. The quiescent halo extends 
across the entire object for ~80 kpc. It is enriched with heavy elements and apparently 
ionized by the continuum from the active nucleus. 

The properties of the quiescent halo and its origin are discussed in this paper. We 
propose that it could be part of a structure that surrounds the entire object, although 
its nature is not clear (a rotating disc? low surface brightness satellites? a cooling flow 
nebula? material ejected in galactic winds? other?). 
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1 INTRODUCTION 

Extended Lya regions are a common feature of high red- 
shift radio galaxies (z >2, HzRG) and quasars (Heckman 
et al. 1991; see also narrow band Lya images of HzRG in, 
e.g, Kurk et al. 2001; Chambers, Miley & van Breugel 1990; 
McCarthy et al. 1990b). Most morphological and kinematic 
studies are based on the high surface brightness regions. 
These regions are clumpy, irregular and often aligned with 
the radio axis. They are characterized by extreme kinema- 
tics, with measured FWHM and velocity shifts >1000 km 
s _1 (Baum & McCarthy 2000, Villar-Martm, Binette & Fos- 
bury 1999, McCarthy, Baum & Spinrad 1996), compared 
to values of ~few hundreds in low redshift radio galaxies 
(Baum, Heckman & van Breugel 1990, Tadhunter, Fosbury 
& Quinn 1989). Although gravitational motions cannot be 
rejected, it is likely that a perturbing mechanism is responsi- 
ble for the extreme kinematics. The apparent connection be- 
tween the radio and the kinematic properties suggests that 
such a mechanism could be shocks generated by the interac- 



tion between the radio outflow and the gas in situ (van Ojik 
et al. 1997). Galactic winds (Heckman, Armus & Miley 1990) 
are an alternative possibility. 

In addition to these regions, low surface brightness Lya 
halos extending beyond the radio structures have been de- 
tected in some HzRG. The kinematic properties of such halos 
have been studied in detail only in one case, the radio galaxy 
MRC1243+036 (van Ojik et al. 1996). The halo shows quies- 
cent kinematics [FWHM(Lya)~ 250 km s _1 ] compared to 
the regions inside the radio structures [FWHM(Lya)~1200 
km s -1 ]. The authors propose that the halo is a large rota- 
ting gaseous disc originating from the accretion associated 
with the formation of the galaxy. 

Such quiescent low surface brightness halos (LSBHs, 
hereafter) are important since they show the gas properties 
unaffected by kinematic perturbations. We are undertaking 
a research program using high S/N Keck spectroscopy of a 
sample of high redshift radio galaxies (z >2.5) whose goal is 
to search for kinematically unperturbed Lya halos in HzRG. 
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We will study the kinematic and ionization properties, as 
well as observed properties such as surface brightness, size 
and luminosity. Constraints will also be set on the halo mass 
and its possible origin. 

The results on the radio galaxy USS0828+193 are pre- 
sented in this paper. In a forthcoming paper we will discuss 
the results on the rest of the sample. 



2 OBSERVATIONS AND DATA ANALYSIS 

USS0828+193 (Rottgering, Miley & Chambers 1995) 
(z =2.57) is a large radio source (~70 kpc)^ showing a dou- 
ble morphology and a jet extending from the core to the 
northern hot spot. The optical HST image shows an irregu- 
lar morphology consisting of several clumps aligned with the 
radio axis (Pentericci et al. 1999, see also Fig. 1 top panel in 
this paper). The high polarization level (~ 10 per cent) im- 
plies that the UV rest frame continuum is dominated by con- 
tinuum from the active nucleus (AGN) scattered by (most 
probably) dust in the extended gas (Vernet et al. 2001). 

The spectra were obtained with the Low Resolution 
Imaging Spectrometer (LRIS, Oke et al. 1995) with its po- 
larimeter (Goodrich, Cohen, Putney 1995) at the Keck II 10 
m telescope in December 1997. We used a 300 line mm -1 
grating and 1 arcsec wide slit which provide a dispersion of 
2.4 A pixel -1 and an effective resolution of FWHM~10.5 
A (instrumental profile). The exposure time was 5 hours. 
The seeing (~1.0 arcsec) was measured with the M5V type 
star located inside the slit. The slit was oriented along the 
radio axis, with PA 44° (Carilli et al. 1997). For a more de- 
tailed description of the observations and data reduction see 
Vernet et al. 2001. 

Each individual frame was calibrated in wavelength and 
corrected for slit curvature using arc spectra. This initial 
wavelength calibration was refined using strong sky lines. In 
order to compensate for residual distortions (shifts) due to 
flexures in LRIS, we extracted from every frame a section 
around Lya+NV 1240 and one around CIV 1550+HeII 1640 
(CIV and Hell respectively, hereafter). The individual sub- 
frames corresponding to the same spectral range were spa- 
tially aligned (using the continuum centroid) and combined. 
The two final frames (one for Lya+NV and another one 
for CIV+Hell) were also spatially aligned, assuming that 
the continuum centroid has the same spatial position at all 
wavelengths. This is reasonable in this small spectral range, 
taking also into account that the spatial resolution element 
has a large physical size (~ 5 kpc). 

The spectra (Lya+NV and Hell+CIV frames) were di- 
vided in several apertures along the spatial direction and a 
spectrum was extracted for each one. The apertures were 
selected so that the gas has apparently similar kinematic 
properties across the spatial extension of a given aperture. 
The apertures were also selected to obtain enough S/N ratio 
in Lya, CIV and Hell (if possible) to fit the line profiles (see 
Figs. 3 and 4). Separate apertures were defined for the gas 
beyond the radio structures. 



* We assume Hq=72 km s 1 Mpc 1 and qo=0.5. In this cosmo- 
logy, 1 arcsec corresponds to ~5 kpc 
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Figure 1. Top: overlay between the optical HST image and the 
radio VLA contours. The 2-dimensional spectrum of the Lya 
emission line (bottom panel) is spatially aligned with the ra- 
dio/optical images. The spatial zero is the position of the con- 
tinuum centroid measured on the Keck spectra. The vertical lines 
indicate the outer edges of the radio lobe. 



3 RESULTS 

3.1 Radio/Optical overlay 

Figure 1 (top panel) shows the overlay between the optical 
WFPC2 HST image and VLA contours (see Pentericci et al. 
1999 for a description). The 2-dimensional Lya Keck spec- 
trum is shown on the bottom panel spatially aligned with 
the image. We identified the radio core with the brightest 
component in the NICMOS HST image. Our main conclu- 
sions would not be affected if we chose to identify the radio 
core with the vertex of the ionization cone described by Pen- 
tericci et al. (1999) or the brightest optical component. 

An interesting feature is already obvious in this figure: 
a low surface brightness Lya halo extending for ~ 16 arcsec 
or 80 kpc inside and beyond the radio structures. It extends 
to much larger distances than the material observed on the 
HST images (~4 arcsec). The HST image shows, therefore, 
only a small fraction of the total distribution of material. 
The kinematics of this halo is more quiescent than the kine- 
matics of the brightest regions. 



3.2 The gas kinematics 

We show in Fig. 2 the 2-dimensional spectrum (contours) 
of the main emission lines. Seven apertures were selected 
in the spatial direction for the kinematic analysis, which are 
indicated in Fig. 3. The corresponding 1-dimensional spectra 
in the Lya and the CIV+Hell regions are shown in Fig. 4. 
Apertures 1, 2, 6, 7 correspond to the low surface bright- 
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Figure 2. Contour plots of the main emission lines. 
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Figure 3. 2-dimensional Lya+NV spectrum and selected aper- 
tures. 



ness Lya halo. Although very noisy, CIV and Hell are de- 
tected in apertures 2 and 6 (see Fig. 4). 

The Lya and Hell spectral profiles were analysed and 
fitted in each aperture with one, two or three Gaussians, 
depending on the quality of the fit. CIV, in spite of being a 
strong line, turned out to be very difficult to constrain due 
to the presence of the doublet components (for a discussion 
on the uncertainties on interpreting the kinematics of the 
extended gas in HzRG using rest frame UV emission lines 
see Villar-Martm et al. 2001). Therefore the results of the 
multiple component fit are only shown for Lya and Hell. 
Those presented for CIV correspond to a single Gaussian 
fit. For each kinematic component, the FWHM (corrected 
for instrumental broadening in quadrature) and the velocity 
shift relative to the emission at the continuum centroid were 
calculated (the Hell line at the continuum centroid was used 
to estimate the zero velocity) . Some of the fits are shown in 
Fig. 5. 

Three kinematic components are found in Lya: (i) a 
narrow component with FWHM<400 km s _1 detected in 
all apertures; (ii) a broad component with FWHM~1200 
km s" 1 detected in all apertures except 6 and 7 (note this 



is also detected beyond the radio lobe in aperture 1); (iii) 
a very broad component with FWHM~2500 km s~ x only 
found in apertures 3 and 4 (brightest regions). The best 
fit to the Hell profile reveals two components in apertures 

3 and 4 similar (in FWHM) to components (i) and (ii) in 
Lya. Narrow Hell (aperture 2) and CIV (apertures 2 and 
6) are also detected. 

(i) The narrow component (<400 km s _1 ): The spatial 
variation of the kinematic properties of this component are 
shown on panels a (FWHM) and b (Vs) in Fig. 6. Diffe- 
rent symbols are used for Lya (circles), Hell (stars) and 
CIV (triangles). This component extends across the whole 
object (~80 kpc) and beyond the radio structures. The kine- 
matic properties are very uniform across the whole spatial 
extension. There are no obvious changes associated with the 
radio structures. 

Lya is very sensitive to absorption by neutral hydrogen 
and the interpretation of the spectral profile in terms of pure 
kinematics is risky. In fact, van Ojik et al. (1997) found 
that the entire blue wing of the Lya (spatially integrated) 
profile of this galaxy is absorbed by neutral gas associated 
with the galaxy. However, the fact that a similar narrow 
component has been found in Hell (which is not a resonance 
line) strongly supports the existence of a narrow component 
in Lya in addition to the broad components. This is also 
supported by the very narrow profile of CIV and Hell in 
apertures 2 and 6 and the very narrow peak observed on 
top of the CIV and Hell lines in the inner apertures, 3 and 

4 (see Fig. 4). 

(ii) The broad component (^1200 km s _1 ) (open circles 
for Lya, open triangles for CIV and stars for Hell in pa- 
nels c and d in Fig. 6). This component is also spatially ex- 
tended (~11 arcsec) and with quite uniform FWHM across 
the object, with no apparent link with the radio structures. 
This kinematically perturbed gas is also found in the Lya 
halo. The measured FWHM is consistent within the errors 
for the three emission lines (Hell, CIV, Lya) where it is de- 
tected, with marginal evidence for narrower Lya. This could 
be due to absorption by neutral hydrogen. This component 
extends beyond the radio structures (aperture 1). Seeing 
effects might be responsible for spreading the detection of 
this component beyond the radio lobes. Astrometry errors 
are not likely to be responsible (see below). 

(iii) The very broad component (^2500 km s _1 ) (solid 
circles in panels c and d in Fig. 6) is spatially extended and 
it is only detected in Lya in apertures 3 and 4. 



4 DISCUSSION 

4.1 The nature of the narrow component 

This component (which extends across the entire object and 
beyond the edge of the radio lobes) shows no apparent asso- 
ciation with the radio structures. The kinematic properties 
are similar to those of low redshift radio galaxies (FWHM 
and velocity shift of a few hundred km s _1 , e.g. Tadhunter, 
Fosbury & Quinn 1989). This suggests that there is no 
(or little) interaction and this gas is unperturbed by the 
jet/radio lobes. Absorption of Lya photons could explain 
the difference in the velocity shift of Lya and Hell in the 
inner apertures. The kinematics of this narrow component 
is consistent with gravitational motions. 
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Figure 4. Spectra extracted for the seven apertures indicated in Fig. 3. The expected position for the CIV and Hell lines has been 
indicated with open arrows in the noisiest spectra. CIV and Hell are detected in the low surface brightness halo (see spectra of apertures 
2 and 6). Therefore, the gas is ionized and contains heavy elements. There is marginal evidence for continuum detection in the halo. The 
interpolation (indicated with 'x') on the red side of CIV and Hell in apertures 1 and 6 was done to remove residuals from sky emission 
lines. 



This gaseous component also emits CIV and Hell lines. 
It is the first time that lines other than Lya are detected 
from the low surface brightness halo of a HzRG. This im- 
plies that the gas is ionized. It is not a 'mirror' of neutral 
hydrogen that reflects Lya by resonant scattering (Villar- 
Martm, Binette & Fosbury 1996). 

There is marginal evidence for continuum detection in 
the halo (see Fig. 4). The nature of this continuum is very 
uncertain. It could be due to scattered light, although stars 
cannot be rejected. The large equivalent widths of CIV and 
Hell suggest that the gas is not ionized by stars, but rather 
by a mechanism related to the active nucleus. This is also 
implied by the CIV/HeII~l value measured for the narrow 
component (at least in aperture 6 and maybe also aper- 
ture 2, see Fig 4), similar to the ratio measured in aper- 
ture 3 and consistent with measurements in other HzRG 



(e.g. Villar-Martm, Tadhunter & Clark 1997). This suggests 
the same ionization mechanism in the low surface bright- 
ness halo and the bright regions, in spite of the different 
kinematic properties. Since the halo extends beyond the ra- 
dio structures where shocks are unlikely to be present, the 
ionization mechanism is probably ionization by the central 
AGN (ionization by the continuum generated by hot shocked 
gas is an alternative mechanism). 

On the other hand, the detection of CIV emission shows 
that the halo has been enriched with heavy elements at very 
large distances from the nuclear region. It is not possible 
to constrain the metallicity using only the CIV/Hell ratio, 
but solar abundances are certainly possible (Villar-Martm, 
Tadhunter & Clark 1997, Vernet et al. 2001). A galactic 
wind, whose presence is suggested by the broad component 
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Figure 5. Examples of some of the fits to the Lya and Hell lines. The fit and the individual components arc shown. 



detected at similar distance from the nucleus, might be the 
mechanism responsible for the enrichment of the gas. 



be scattered light from a hidden broad line region (Vernet 
et al. 2001). 



4.2 The nature of the broad (ii) and very broad 
(iii) components 

The broad component (~1200 km s _1 ) is apparently de- 
tected beyond the Eastern radio lobe. Errors in the Astro- 
metry are unlikely to be responsible: if the active nucleus 
were located at the vertex of the ionization cone described 
by Pentericci et al. (1999) or at the brightest optical feature, 
the broad component would extend even further beyond the 
edge of the radio lobe. However, we do not discard the possi- 
bility that seeing effects might be smearing the emission 
enough that the broad component actually does not extend 
beyond the radio lobes. 

If the detection beyond the radio structures is real, this 
implies that the motions of this gas are not a consequence 
of jet-cloud interactions. Evidence for perturbed gas beyond 
radio structures has been found in other objects (e.g. Villar- 
Martm, Binette and Fosbury 1999, Solorzano-Iharrea, Tad- 
hunter & Axon 2001). Galactic winds could be an alternative 
possibility. Clouds exposed to ram pressure of the wind can 
be accelerated to several hundred km s _1 and shifted by 
~1000 km s" 1 (Heckman, Annus & Miley 1990). The inte- 
gration along the line of sight through several such clouds 
could explain the broad FWHM of this component. Gas free- 
falling into the galaxy during the formation process from a 
large distance can also explain FWHM of this order (Lehn- 
ert & Becker 1998). The nature of the very broad component 
(~2500 km s" 1 ) in apertures 3 and 4 is not clear. It might 



4.3 Mass and density of the ionized gas in 
USS0828+193 

(Hereafter, we will refer to the gas emitting the narrow com- 
ponent as the low surface brightness halo, LSBH). 

The integrated Lya flux Fi ya of the LSBH along the 
slit is 1.64xl0 -15 erg s _1 cm -2 , which corresponds to a 
luminosity L_L ya =3.8xl0 43 erg s _1 . For pure case B recom- 
bination and T =10 4 K, L-Lyu is given by (McCarthy et al. 
1990a): 



-'Lya 



= 4 x 10 



-24 2 c T/ 

n e f V erg s 



(1) 



where / is the volume filling factor of the gas, n e the 
electronic density and V is the total volume inside the slit. 

We have calculated the volume assuming that the gas 
of the LSBH is ionized by the active nucleus (see §4.1) and is 
homogeneously distributed. In this scenario, the ionized halo 
has a biconical geometry. We have assumed that the cone 
axis is on the plane of the sky with an opening angle for the 
individual cones of ~90° (Barthel 1989). The 2-dimensional 
spectra show an extension of ~8 arcsec (40 kpc) along the 
slit from the AGN. The volume of the bi-cone slice included 
in the slit (1 arcsec wide or 5 kpc) is ~3.2xl0 68 cm 3 . Assu- 
ming / =10~ 5 (McCarthy et al. 1990a), we obtain n e ~50 
cm~ 3 . This is a very crude estimate. Large uncertainties 
are due to: a) possible scatter on the / values. Values for 
/ of the order of 10 _7 -10 -8 have been estimated for Lya 
nebulosities associated with high redshift quasars (Heckman 
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Figure 6. Kinematic properties of the individual components de- 
tected in the different apertures. The FWHM appears on the left 
panels and the velocity shift relative to the Hell emission at the 
continuum centroid is plotted on the right panels. The narrow 
component is plotted in panels a and b and the broad and very 
broad components are plotted in panels c and d. The dashed ver- 
tical lines mark the outer edge of the radio lobes. Arrows indicate 
upper limits. The horizontal lines in panel a show the spatial 
extension of the aperture indicated with a number. Notice the 
uniform properties of the narrow component across the whole ob- 
ject and the similarity with low redshift radio galaxies kinematic 
properties. The detection of the broad component (panels c and 
d) beyond the radio structures is also interesting, since it suggests 
that jet-cloud interactions are not responsible for the perturbed 
motions. 



et al. 1991) and some low redshift radio galaxies (Clark et 
al. 1998); b) the fact that the measured Lya luminosity is 
probably lower than the case B value (due to the dust and/or 
neutral hydrogen absorption) and c) the assumptions on the 
geometry. The estimated n e value is likely to be a lower 
limit for the density and it could be at least one order of 
magnitude higher. 

The mass of ionized gas in the LSBH is given by 



M{H + ) = n e xV xm p x f 



(2) 



where V is in this case the total volume of ionized gas 
and m p the mass of the proton. (We assume that the den- 
sity is constant across the halo for simplicity, but it is possi- 
ble that the density increases inwards). Considering the two 
cones we obtain a total volume of 2.3xl0 69 cm 3 and there- 
fore M(H + ) ~9.6xlO s M a , for n e =50 cm~ 3 . This is proba- 
bly un upper limit. Since M(H + ) oc f 1 ^ 2 and taking into 
account the uncertainties on /, the mass could be at least 
ten times lower. 

The halo material outside the ionization cones is likely 



to be neutral (or have a low ionization level due to, for 
instance, metagalactic radiation). Evidence for reservoirs of 
neutral gas associated with distant radio galaxies has been 
found by van Ojik et al. (1997). The authors found that the 
Lya profile in eleven (including USS0828+193) out of eigh- 
teen distant radio galaxies shows deep troughs which they 
interpret as HI absorption by absorbers that are likely to 
be physically associated with the galaxy hosting the radio 
source or its direct environment []. 

The LSBH in USS0828+193 (and other HzRG) could be 
part of a gaseous reservoir that surrounds the object com- 
pletely. The gas outside the ionization cones is neutral and 
we see it in absorption and the gas inside the cones is ioni- 
zed and we see it in emission. If we assume that the neutral 
gas of the halo has the same density and filling factor as 
the ionized fraction, we obtain M(H°) ~2.4xl0 9 M Q . The 
mass of ionized (T ~10 4 K) and neutral gas in the halo is 
therefore M H++H o ~3.4xl0 9 M Q . 

Following a similar procedure, we obtain n e ~ 150 cm -3 
for the broad line emitting gas (the same / and cone open- 
ing angle were used). To calculate the volume, two cones 
of height 4 and 6 arcsec respectively were considered. These 
are the projected distance values at which the broad compo- 
nent is detected at both sides of the continuum centroid (see 
Fig. 6). This gives a volume of ~1.3xl0 68 cm 3 inside the slit 
and 1.2xl0 69 cm 3 as the total volume inside the cones. The 
mass of broad emission line gas inside the ionization cones 
is then ~1.5xl0 9 M . 

We conclude that a large fraction of the gaseous (with 
T <10 4 K) mass in HzRG might be in large halos that su- 
rround the entire object. Part of the halo (the gas inside 
the ionization cones) is ionized and we see it as low surface 
brightness line emission. The rest is neutral (or with very 
low ionization level) and we see it as absorption features 
imprinted on the Lya profile 



4.4 The origin of the LSBH in USS0828+193 

We have considered several scenarios: 

• A rotating disc (evidence for large discs associated with 
low redshift radio galaxies is discussed in §4.5.1). Signs for 
rotation of the extended gas in nearby radio galaxies have 
been found by several authors (Heckman et al. 1985, Baum, 
Heckman & van Breugel 1990, Tadhunter, Fosbury & Quinn 
1989). In this case, M rot = G R S ^ p where R is the radius 
of the disc (40 kpc), V is half the amplitude of the rotation 
curve and i is the inclination angle of the disc with respect 
to the plane of the sky. This formula is applicable when we 
measure the rotation along the line of nodes. Velocity curves 
that resemble rotation have been found in several HzRG, but 
it is not clear whether this is an effect of multiple kinematic 
components (Villar-Martfn et al. 2001). 

The current data do not allow us to set any constraint on 



t An interesting case is the radio galaxy MRC0943-242 
(z =2.92), which shows large scale absorption troughs (Rottgering 
et al. 1995, Binette et al. 2000), both in Lya and CIV. Binette 
et al. concluded that the absorbing gas is metal poor (Z ~ 0.01 
Zq) and is located much further outside the zone of influence of 
the radio jet cocoon 
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Figure 7. Same as panel b in Fig. 6 but with the measured ve- 
locity shift for the very narrow peak observed in CIV in the inner 
apertures (solid triangles). CIV follows a similar pattern as Lya. 



the mass if the halo is a rotating disc: we do not have infor- 
mation on the inclination angle of the disc (imaging would 
help to constrain it), neither on the direction of the line of 
nodes, which is likely to be misaligned with respect to the 
slit position. Van Ojik et al. (1996) propose that the giant 
Lya halo associated with the radio galaxy MRC 1243+046 
(z =3.6) has settled as a rotating disc whose orientation (the 
plane of the disc) is within ~20° of the orientation of the 
radio axis. However, this is contrary to the trend observed in 
low redshift powerful radio galaxies with signs of rotation. 
In such objects, the rotation axis of the gas and the radio 
axis trend to be aligned within 20° (Heckman et al. 1985) . If 
this is the case for the LSBH in USS0828+193, the rotation 
axis would be approximately along the slit, since this was 
aligned with the radio axis. 

The Lya velocity curve of the narrow component in 
USS0828+193 does not resemble a standard Keplerian rota- 
tion curve for a slit located along the line of nodes. However, 
it is interesting to note that the observed pattern is similar 
to what we expect from a rotating disc if the slit was not 
aligned with the line of nodes, but rather close to the per- 
pendicular direction and with some impact parameter with 
respect to the rotation center (i.e. there is a distance be- 
tween the center of the slit and the center of rotation). Such 
a pattern is expected if the disc has a certain inclination rel- 
ative to the plane of the sky. A close analysis of the CIV line 
profile in the two inner apertures (3 and 4, see Fig. 4) reveals 
the presence of a very narrow peak (on top of a dominant 
broader component). The velocity shift of this narrow peak 
has been plotted in Fig. 7 (solid triangles); CIV shows a 
pattern similar to Lya. 

3-dimensional spectroscopy would be very valuable to ob- 
tain information about the geometry and the kinematic pat- 
tern along different directions across the halo. 

• A virialized system: The halo could be a virialized sys- 
tem that consists of N individual clouds of mass m that 



move with random motions 0. This could represent the sce- 
nario where the halo is formed by galactic satellites. Ap- 

5 R V, 2 



plying the virial theorem Md yn 



(Carroll & Ostlie 



1999), where R is the radius of the system (40 kpc) and V v 
is the radial velocity dispersion of the clouds within the halo 
(given by the a of the emission lines ~130 km s -1 ) we obtain 
Mdyn = 8X10 11 Mq (this is the total mass of the system, 
luminous+dark) . 

• Inflow: In this case M in j = v £ , where Vr is the in- 
falling velocity of a particle at distance R from the cen- 
tral mass Mi„f. We have assumed that the axis of the ioni- 
zed cones (opening angle 0=90°) is on the plane of the sky. 
The emission at each spatial position along the slit results 
from the integration of the gas emission within the cone 
slice inside the slit and along the line of sight (see Fig. 8). 
In such case, the FWZI of the narrow emission line at a 
given slit position gives the maximum projected velocity dis- 
placement between the gas at the two opposite edges of the 
slice cone in that spatial position, i.e., FWZI = V02 — V01 = 
2x V R sin{0/2), where V&a = -V01 = Vr sin(0/2). On the 



other hand, R = - 
radius. Therefore 



WT2) 



, where r is the observed projected 



M, 



inf 



r FWZI a 1 

cos(6/2) ( 2 sin(6/2)' 2 G 



(3) 



In USS0828+193, FWZI~430 km s" 1 at r =40 kpc (8 arc- 
sec) as measured from the spectrum. We obtain M in f ~ 
6xlO n Mq for the total mass of the system. 

Values of this order have been estimated for the masses 
of low redshift radio galaxies and are smaller than typical 
masses of cD elliptical galaxies at low redshift (Tadhunter, 
Fosbury & Quinn 1989). If this result is confirmed, we would 
conclude that USS0828+193 will not become like the cD 
galaxies we see at the present epoch, unless mergers are in- 
volved. 

Is it possible that the material in the halo is falling to- 
wards the center in a cooling flow manner (rather than in 

free fall) ? Q. We have investigated whether this scenario is 
plausible. We have assumed that the emission line gas is in 
pressure equilibrium with the hot phase. Typical tempera- 
tures of Tdouds =10 000 K and That =10 r for the optical 
emitting gas and the hot phase respectively and n c i ou d s =50 
cm -3 for the emission line gas (§4.2), give rih ot =0.05 cm -3 
(the uncertainties on n c i ou d a and Thot allow only crude esti- 
mates in the discussion that follows). 

Adopting a value for the cooling function 7(T=10 7 K) 
= 5xl0~ 23 erg cm 3 s _1 (Raymond, Cox & Smith 1976), 
and combining the three equations in §3 (Nulsen, Stuart 
& Fabian 1984) the mass of hot gas in the halo (inside the 
visible radius of 40 kpc) of USS0828+193 is ~4x 10 s Mq and 
its radiative cooling time is ~6x 10 7 years. This cooling time 
is rather short (much shorter than the age of the Universe 
at that time). On this basis, it is possible that the gas we 
observe from the LSBH has cooled from the hot phase. 



t Pettini et al. (2001) also calculate dynamical masses of Ly break 
galaxies in terms of gas clouds. 

§ In the free fall scenario, the only force acting on the gas is 
gravity. In the cooling flow scenario, the pressure gradient in a 
hot corona acts against gravity 
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Calculating the total mass in the infalling scenario 
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Figure 8. Infall scenario: Vo\ and V02 are the velocities of par- 
ticles in positions 1 and 2, projected on the line of sight (l.o.s) 
and Vr is the radial (infalling) velocity towards the center. 



The mass that has cooled down is ~10 9 Mq (considering 
only the ionized gas in the LSBH [see §4.2]) . If this happened 
m ~6xl0 7 yr, the mass deposition rate is 17 Mq yr 1 (va- 
lues of this order have been measured in nearby cooling flow 
galaxies, Heckman, Armus & Miley 1990) 

We have assumed that the density of the halo is constant. 
However, this is a very simplified scenario and the density 
profile is likely to be centrally condensed (Haiman and Rees 
2001, HR2001 hereafter). The cooling time varies as n" 1 
(Nulsen, Stewart & Fabian 1984), therefore, the outer, less 
dense gas will need longer time to cool down. As van Ojik et 
al. found, when we account for this effect, the cooling time is 
still rather short at large distances (r=100 kpc) for a density 
behaviour of both oc r" 1 and oc r" 2 (2xl0 8 and 2.5x10 s 
years respectively). There has been, therefore, enough time 
for the gas to cool down. 

In a cooling flow scenario, we expect very narrow emission 
lines at large distances (<100 km s _1 at >tens kpcs, Fabian 

et al. 1987) a chaotic pattern in the velocity curve and 
a trend of the line width to decrease with increasing radius 
(Heckman et al. 1989) . This behaviour is not obvious on the 
data but we are limited by the low spectral resolution. The 
kinematic properties of the halo do not help to constrain the 
validity of the cooling flow scenario. Higher resolution would 
help to measure the FWHM more accurately at different 
spatial position. 

• Outflow: the halo could have been deposited by massi- 



^ Emission lines as broad as several hundred km s _1 are observed 
in low redshift cooling flow nebulae (e.g. Fabian et al. 1987, Heck- 
man et al. 1989). The authors propose that galaxy interactions 
may generate turbulence. 



ve outflows (winds could explain the existence of heavy ele- 
ments in the halo). We have calculated some basic param- 
eters characterizing a superwind which could explain the 
kinematic properties of the LSBH in USS0828+193. The 
energy injection rate for an energy-conserving bubble in- 
flated by energy injected at a constant rate and expanding 
into a uniform medium with an ambient density no (cm -3 ) 
is related to the radius (in kpc) and velocity (in 10 2 km s" 1 ) 
by (Heckman, Armus & Miley 1990) 



dE/dt = 3x 10 41 r 2 pc Wi 00 no ergs 1 



(4) 



and the dynamical time scale of the nebula is given by 



tdynfi ~ 10 r kpc /v 100 Myr 



(5) 



where V100 would be halve the separation in velocity be- 
tween the two components of the double peaked profile pro- 
duced by the expansion. 

The resolution of the spectrum is too low to resolve the 
double peak that the expansion of the bubble would pro- 
duce. We have estimated that two unresolved emission lines 
separated by a given velocity 8v would produce a line of 
FWHM > 2x5v — A xuioo (all velocities in units of 100 km 
s" 1 ). Therefore, from the measurements in our spectrum 
(a ~130 km s" 1 ) we obtain uioo <0.76. 

The estimated density n e =50 cm -3 is probably the den- 
sity of the post-shock gas in the expanding bubble, no in the 
dE/dt equation is, however, the pre-shock density, which is 
likely to be much smaller (^several cm -3 , Heckman, Ar- 
mus & Miley 1990). For n e =l cm" 3 the implied values of 
tdyn and dE/dt are ~530 Myr (the corresponding time scale 
for ARP 220 is ~100 Myr, Heckman, Armus & Miley 1990) 
and 2.1~10 44 erg s _1 respectively. The required star forming 
rate is SFR~500 Mq and a supernova rate of ~4 SN yr -1 . 
Typical SFR values for Ly break galaxies (dust corrected) 
are several hundred Mq yr -1 (e.g. Sawicky & Yee 1998). 

The expected total infrared luminosity of such a system 
would be Lir ~3xl0 12 L . This would place USS0828+193 
in the category of ultraluminous infrared galaxies. 

The kinematic signature of an expanding bubble (or bi- 
conical structure) would be a region of double peaked emis- 
sion line profiles. Higher spectral resolution might reveal this 
pattern. 



4.5 Is this a halo from which the galaxy is still 
forming ? 

As Haiman & Rees discuss (2001), in the process of forma- 
tion of a galaxy, smooth gas collapses from the virial ra- 
dius to its final orbit radius. Numerical simulations have 
revealed a more complex process, where a fraction of the in- 
falling gas forms smaller clumps early on; these clumps then 
progressively merge together, collide and dissipate to form 
larger systems. In these scenarios, a robust feature of galaxy 
formation is at least in the early stages a spatially extended 
distribution of gas. HR2001 showed that if the halo is illumi- 
nated by ionizing radiation from a central quasar, it should 
be observed as an extended Lya fuzz with surface brightness 
~10" 17 erg s" 1 cm" 2 arcsec -2 . 

We compare in this section the properties of the LSBH 
in USS0828+193 with those expected for the extended gas 
condensation in the early stages of galaxy formation models. 
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• Total mass: At z =3 (total) halo masses in the range 
~4xl0 10 -10 13 Mq are expected for a virial temperature in 
the range T vir =2 x 10 5 K - 2 x 10 7 K (HR2001). We have 
obtained Mtot ^-SxlO 11 Mq for the two scenarios (infall 
and virialized system of clouds) where we have been able to 
calculate the total mass of the system. This value lies in the 
range predicted by the models. 

• F(Lya)ss: The expected surface brightness for the 
Lya fuzz is in the range ~10 -18 -10 -16 erg s _1 cm" 2 
arcsec -2 (HR2001). This is consistent with the measured 
value in USS0828+193 ~ 3xl0~ 17 erg s" 1 cm -2 arcsec" 2 . 

• Linear size: The linear spatial extension of the Lya 
emitting gas would be a fraction of the virial radius R V i r ~ 
10 - 100 kpc (HR2001), consistent with the linear scale of 
the halo (radius) in USS0828+193, ~40 kpc. 

• Black hole mass: According to HR2001 model, the 
minimum mass required to ionize the halo is given by: 

M 6h ~6xlO 8 M 0(T J^-f/ 3 (i±i) 4 (6) 

Assuming M to t=8xl0 11 Mq as above, we find a lower 
limit for the black hole mass of ~5x10 7 Mq. This value is 
somewhat lower than the black hole mass we expect accor- 
ding to Gu, Cao and Jiang (2001), who found that the vast 
majority of the quasars in the Uy, S4 and S5 catalogues 
have typical black hole masses of >10 8 Mq. 

• Heavy elements: One of our most important results 
is the finding of heavy elements in the LSBH. Any model of 
protogalactic halos should consider metal enrichment. This 
is mentioned by HR2001, although not discussed in details. 

The properties of the low surface brightness halo in 
USS0828+193 are, therefore, consistent with the expected 
properties for the extended gaseous halo, a required ingre- 
dient for galaxy formation models. 

However, the existence of the halo does not necessarily 
imply that USS0828+193 is a galaxy in the process of forma- 
tion. Similar halos have been found in lower redshift radio 
galaxies, where the galaxy is expected to be fully formed. 
An example is the radio galaxy 3C368 at z —1 (Stock- 
ton, Ridgway & Kellogg 1996). This object shows low sur- 
face brightness [Oil] 3727 emission with apparently quies- 
cent kinematics, in addition to the bright kinematically per- 
turbed structures. Other examples at lower redshift are 3C34 
(2=0.69), 3C435A (2=0.47) (see 2-dimensional [Oil] 3727 
spectra in Solorzano-Inarrea, Tadhunter & Axon (2001). The 
halos extend for more than 100 kpc). Such halos might be 
similar to the LSBH we have discovered in the Lya light 
in our more distant radio galaxies. If this is the case, the 
existence of LSBH associated with distant radio galaxies 
does not imply an early stage in the formation process. On 
the other hand, the consistency of the LSBH properties in 
USS0828+193 with expectations for the initial gas reservoir, 
suggests that the halo gas is part of the material from which 
the galaxy started to grow. 

4-5.1 Comparison with other objects 

We discuss here several structures found around other types 
of galaxies at different redshifts which might have the same 
origin as the low surface brightness halos found around dis- 
tant radio galaxies. 



• Galactic envelopes: observations of absorption line sys- 
tems in the spectra of background quasars have provided 
evidence for large (R ~100 h~ x kpc) extended gaseous en- 
velopes that surround galaxies of a wide range of luminosity 
and morphological type (e.g. Chen, Lanzetta & Webb 2001; 
Lanzetta et al. 1995). Chen, Lanzetta & Webb (2001) have 
proposed accreting satellites as the origin of these structures 
to explain the chemical enrichment at large distances. 

• Giant Lya halos associated with Ly break galaxies: Stei- 
del et al. (2000) found two giant (physical extent >100 
diffuse Lya emitters apparently associated with previously 
known Ly break galaxies at < z >=3.09. The origin of these 
structures is not clear. Proposed explanations are a cooling 
flow nebula (Steidel et al. 2000) or 'hyperwinds' (Taniguchi, 
Shioya & Kakazu, 2001). 

• Large scale (>several tens of kpcs) HI disc like struc- 
tures have been found in several low redshift radio (Morganti 
et al. 2002) and elliptical galaxies (Oosterloo et al. 2002). 
The authors propose that these structures are the results of 
mergers. The LSBH found around USS 0828+193 and other 
distant radio galaxies, might be the progenitors of the discs 
found at low redshift. 

Therefore, giant gaseous halos (> several tens of kpc) 
are found to be often associated with galaxies in general: 
active and non active, of different morphological types and 
luminosities and at different epochs. The origin and nature 
of such halos is not clear and different explanations have 
been proposed for different objects, but the definitive answer 
is uncertain. An interesting possibility is that these halos 
have a similar origin related to the formation process of the 
galaxies. 



5 SUMMARY AND CONCLUSIONS 

By means of the spectral decomposition of the emission line 
profiles in USS0828+193, we have isolated the emission from 
a giant (~80 kpc) reservoir of apparently kinematically un- 
perturbed gas, whose kinematics is consistent with being 
gravitational in nature. Emission lines other than Lya (CIV, 
Hell) have been detected for the first time in such a halo. 
This implies that the gas is ionized (probably by the con- 
tinuum from the AGN) and enriched with heavy elements. 
We find marginal evidence for continuum detection from the 
halo. 

We propose that the low surface brightness halo in 
USS0828+193 (and other HzRG) could be part of a gaseous 
reservoir that surrounds the entire object. We have consi- 
dered several possible scenarios to explain the origin of the 
halo. The gas could have settled in a rotating disc. This is 
suggested by the velocity pattern of the halo emission, al- 
though 3-dimensional spectroscopy would be necessary to 
set tighter constrains on the kinematic patterns. Evidence 
for giant gaseous discs has been found by other authors in 
low redshift radio galaxies and elliptical galaxies. 

The halo could also be a group of individual clouds 
(galactic satellites?) in a virialized system. In this case, the 
expected mass is ~8xlO n Mq. A similar value is estimated 
if the halo is infalling towards the center of the potential 
well. Such values are similar to estimates of masses of low 
redshift radio galaxies, and smaller than typical masses of 
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cD elliptical galaxies at low redshift. If this result is con- 
firmed, we conclude that USS0828+193 will not become like 
the cD galaxies we see at the present epoch, unless mergers 
are involved. 

Alternative scenarios that we have discussed (also con- 
sistent with the observations) are a cooling flow from a hot 
phase and outflows (galactic winds). In this last case, star 
forming rates ~500 M© yr _1 are required, consistent with 
estimates for Lyman break galaxies. 

The properties of the low surface brightness Lya halo 
detected in the radio galaxy USS0828+193 are consistent 
with the expectations for the original gaseous reservoir from 
which the galaxy started to form. This does not necessarily 
imply that the galaxy is in the process of formation. 

We have also discussed the possible link of the low sur- 
face brightness halos with a) the very extended envelopes 
found around galaxies of different morphological types and 
luminosities based on studies of absorption line systems in 
the spectra of background quasars b) the giant Lya halos 
associated with some Ly break galaxies c) the very extended 
reservoirs (claimed to be discs) of neutral gas associated with 
some nearby radio galaxies and elliptical galaxies. 
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